PASJ: Publ. Astron. Soc. Japan , 1-??, 2002 October 25 
© 2008. Astronomical Society of Japan. 



3-D Kinematics of Water Masers in the W51A Region 
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Abstract 

We report proper motion measurements of water masers in the massive-star forming region W 51A and 
the analyses of the 3-D kinematics of the masers in three maser clusters of W51A (W51 North, Main, and 
South). In W 51 North, we found a clear expanding flow that has an expansion velocity of ~70 km s _1 
and indicates deceleration. The originating point of the flow coincides within 0".l with a silicon-monoxide 
maser source near the HII region W 51d. In W51 Main, no systematic motion was found in the whole 
velocity range (158 km s _1 < Vlsr < —58 km s -1 ) although a stream motion was reported previously 
in a limited range of the Doppler velocity (54 km s _1 < Vlsr < 68 km s -1 ). Multiple driving sources of 
outflows are thought to explain the kinematics of W51 Main. In W51 South, an expansion motion like a 
bipolar flow was marginally visible. Analyses based on diagonalization of the variance-covariance matrix 
of maser velocity vectors demonstrate that the maser kinematics in W 51 North and Main are significantly 
tri-axially asymmetric. We estimated a distance to W51 North to be 6.1±1.3 kpc on the basis of the model 
fitting method adopting a radially expanding flow. 

Key words: masers — starsdormation — stars:distances — stars individual (W51A, W51 North, 
W 51 Main, W51 South) — ISM:jets and outflows 



1. Introduction 

Water maser emission is one of the most important 
phenomena in the study of star formation, often based 
on data obtained using very long baseline interferom- 
ctry (VLBI) with high angular and velocity resolution 
(e.g., ©00 [cite]cite.rei81Rcid & Moran(1981); ©00 
[cite]cite.eli92Elitzur(1992)). Analyses of spatial po- 
sitions, Doppler velocities, and proper motions of 
individual maser features with a typical size of 1 
AU (©00 [cite]cite.rei81Reid & Moran(1981)) have re- 
vealed the 3-D gas kinematics around young stellar 
objects (YSOs) (e.g., ©00 [cite] cite. gen81aGenzel et al. 
(1981)a, hereafter G81; ©00 [cite]cite.gen81bGenzel 
|00 [cite]cite.sch81Schncps ct al. 

S81; 030 [cite]cite.rei88Reid 
[cite]cite.gwi92Gwinn, Moran, 
|00 [cite]cite.cla96Claussen et al. 
citc.furOOFuruya et al. (2000); ©00 
et al. (2000), hereafter Paper I; 



ct al. (1981)b; 
(1981), hereafter 
et al. (1988); 
& Rcid(1992); 
(1996); §)0 [cite 
[cite! citc.imaOOImai 



30 [cite] cite. torOlTorrcllcs ct al. (2001)). In practice, 
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details of the gas kinematics are complicated but seem 
to depend mainly on the evolutionary status of YSOs. 
Water maser sources sometimes enable measurement 
of the internal motions of giant molecular clouds by 
measuring relative bulk motions between clusters of 
water masers that are very close to each other (Paper 
I; ©00 [cite]citc.tor01Toi-rellcs ct al. (2001)). Such bulk 
motions may be owe to, e.g., propagation of shock layers 
from newly-formed HII regions, cloud contraction by the 
self-gravitation due to the huge mass of a giant molecular 
cloud. 

The massive-star forming region W 51A contains 
at least five independent clusters of water masers: 
W51 North, West, Main, and two clusters in South 
(e4 and c3). All of the maser clusters are inde- 
pendently associated with HII regions and dense 
and young cloud cores exhibiting several species 
of molecular emission (©00 [cite]cite.gen78Genzel 
et al. (1978); 030 [cite]cite.dow79Downes et al. 
(1979); G81; S81; |o|DO [cite]cite.gau87Gaume & 
Mutel (1987); ©00 [cite]cite.gau93Gaume, Johnston, 
& Wilson(1993); 000 [cite]cite.zha95Zhang & Ho(1995), 
1997; ©00 [cite]cite.zha98Zhang, Ho, & Ohashi(1998); ©00 
[cite] cite. lep98Leppanen, Liljestrom & Diamond(1998), 
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hereafter LLD; ©00 [cite]cite.eis02Eisner et al. (2002), 
hereafter EGHMM). In W51N, silicon-monoxide masers 
have been detected ([J)0 [cite] cite. mor92Morita et al. 
(1992); EGHMM). Proper motions of these water masers 
were measured two decades ago (G81; S81). However, 
details of the kinematics are still obscured because of 
limited numbers of measured proper motions (<30, S81; 
G81). In W51 North, a separation motion has been 
recognized between t wo d ominant " sub-clusters" of maser 
features (see Sect. BT), so-called the "NW Cluster" 
and the "Dominant Center Reference Cluster" (S81). 
EGHMM reconfirmed the similar separation motion in 
the patterns of the sub-clusters, which are similar to 
bow shocks, by comparing the patterns observed in 1998 
with those of S81. In W51 Main, LLD found a stream 
from the "Double-Knot" sub-cluster in the south-east 
direction. The observed velocity range, however, was very 
limited (54 km s _1 < Ulsr < 68 km s™ 1 ). Throughout 
the entire range (—50 km s" 1 < Vlsr < 130 km s _1 ) the 
kinematics are predominated by random motions (G81). 
In other maser clusters, their kinematics have never been 
well understood. 

Here, we report monitoring observations of the W 51A 
water masers with the Japanese domestic VLBI network 
(J-Net) [] (§00 [cite]cite.omo940modaka et al. (1994)). 
Sect. 2 describes the VLBI observation and data reduc- 
tion. Sect. 3 summarizes the revealed 3-D kinematics of 
the individual clusters of water masers. Sect. 4 discusses 
the origins and related issues of the maser kinematics. The 
estimation of the distance to W 51A and relative 3-D bulk 
motions among the maser clusters are also described. 

2. Observations and Data Reduction 

The VLBI observations were made at five epochs for 
a period of 8 months in 1999, using three or four J-Net 
telescopes. Table |l| gives the status of the J-Net obser- 
vations. The spatial resolution of our observations was 
typically 3 milliarcseconds (mas). W51 North/West, 
Main, and two regions of W51 South, close to the HII 
regions W51 e4 or e8 and e3 (©00 [cite]cite.gau93Gaume, 
Johnston, & Wilson(1993); ©00 [cite]cite.zha98Zhang, Ho, 
& Ohashi(1998)), are located within a 70"-field. We ob- 
served simultaneously the four clusters of masers within a 
single antenna beam (the minimum size of 72" [FWHM] 
at 22 GHz in the array telescopes). At each epoch, an 
observation was made for 10-11 hours including scans of 
calibrators (NRAO 530 and JVAS 2145+067. The re- 
ceived data were recorded with the VSOP terminal (©00 
[cite] cite. kaw94Kawaguchi et al. (1994)) in a single base 
band channel with a band width of 16 MHz, which corre- 
sponds to a velocity coverage of 216 km s" 1 (158 km s _1 
< 14sr < -58 km s^ 1 ). 

The data correlation was made with the Mitaka FX 



1 J-Net includes the 45-m telescope of Nobeyama Radio 
Observatory (NRO), which is a branch of the National 
Astronomical Observatory, an interuniversity research institute 
operated by the Ministry of Education, Culture, Sports, Science 
and Technology. 



correlator (§00 [cite]cite.chi91Chikada et al. (1991)). 
Correlations were performed four times by shifting the 
phase-tracking centers to the locations of the above four 
clusters of masers. An average correlation time was set to 
0.5 s for W51 North and 1 s for W51 Main and two regions 
of W 51 South. The correlation outputs consisted of 1024 
velocity channels with a velocity spacing of 0.21 km s _1 
each. 

Data reduction was made with the NRAO 
AIPS package using normal procedures (e.g., ©00 
[cite]cite.dia95Diamond(1995)). Fringe fitting and self- 
calibration procedures were performed for a Doppler 
velocity channel including a bright maser spot (velocity 
component) in W51 North, also given in Table [|. The 
solutions were applied to all of the data, then maser 
maps were made for all of the maser clusters. The 
typical size of the synthesized beam was 3 mas in the five 
observations (see Table 0). The relative position accuracy 
of a maser spot ranged over 0.02-0.6 mas depending on 
a signal-to-noise ratio and spatial structure of the spot. 
Identification of a water maser feature was made in the 
same procedure shown in several previous papers (e.g., 
Paper I). A relative position accuracy of a maser feature 
was ranged over 0.05-0.6 mas. Relative proper motions 
were measured for maser features identified in at least 
two of the five epochs. 

3. Results 

3.1. Summary of proper motion measurements 

Figure [l] shows several examples of measured relative 
proper motions of water maser features. Maser features 
fundamentally seem to move with constant velocities. The 
deviations from fit lines assuming constant velocity mo- 
tions are within several tenths of a milliarcsccond. Some 
features have a large deviation from the fit lines because 
they are located together with other features within a 
small range, 1 mas in space and 1 km s _1 in velocity, in 
which we were not able to correctly trace the same feature 
from one epoch to another. 

Tables |, |§ and | give parameters of maser features 
with measured proper motions. The numbers of measured 
proper motions were 123, 48, and 10 in W51 North, Main, 
and South (c4), respectively, which are larger than those of 
previous observations for W51 North by S81 and for W51 
Main and South (e4) by G81. An important difference 
between the previous and the present measurements is 
the difference in time separations between the successive 
observing epochs: two years and, at minimum, only one 
month, respectively. 

Even with a much shorter time separation of our obser- 
vations, it has not been possible to measure proper mo- 
tions in a large fraction of detected maser features (> 50 
%) mainly because of growth and decay of maser features 
among the observing epochs. These results imply that 
many of the maser features have lifetimes shorter than 
1-2 months. Unfortunately, no maser proper motion has 
been identified in the W 51 South (e3) region. 

Usually, each of water maser clusters consists of several 
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groups of maser features with a size of 100-1000 AU. In 
this paper, we define such a group of features as a "sub- 
cluster" . 

3.2. W51 North 

3.2.1. Overview 

Figure presents the angular distribution of water 
masers in W51 North and West, which revives those found 
by previous observations (000 [cite]cite.gen78Genzel et al. 
(1978); 030 [cite]cite.dow79Downes et al. (1979); S81; 
EGHMM) . Figure presents the detailed 3-D motions of 
water masers around the two dominant sub-clusters, the 
red-shifted (north-west) and the blue-shifted (south-east, 
Dominant Center Reference) sub-clusters, which clearly 
exhibit a bipolar expanding flow. We found a "bow 
shock" pattern in the SE sub-cluster, which opens toward 
the NW sub-cluster and has been found since 1977 (S81; 
EGHMM). On the basis of the locations of the two sub- 
clusters and this bow-shock pattern, we estimated a loca- 
tion of SiO maser emission in W51N on our map with an 
uncertainty of less than 100 mas, whose position was mea- 
sured by EGHMM with respect to the water masers with 
an uncertainty of less than 50 mas. The originating point 
of the outflow seems to be located around the middle of the 
two maser sub-clusters and to be roughly coincident with 
the location of the SiO maser. While only a separation 
motion between the two sub-clusters has been confirmed 
by S81 and EGHMM, the present result reveals that the 
individual sub-clusters themselves are also expanding. 

On the other hand, motions of maser features far from 
the two sub-clusters do not exhibit any systematic mo- 
tion but random motions with velocities up to 100 km s _1 . 
Some peculiar motions arc found in the two dominant sub- 
clusters and are not due to misidentification of the proper 
motions. Thus, although the expanding flow has been 
visible, the kinematics of the W51 North region is heavily 
disordered dynamically by fast random motions. The de- 
tails of the maser kinematics around W51 West were also 
obscured because only two proper motions were measured. 
It is difficult to compare the maser distribution with that 
of S81 because of too sparse time separation (~20 yr). 

3.2.2. Model fitting the maser kinematics 

In order to estimate kinematical parameters of the ex- 
panding flow and a distance to W51N, we made model- 
fitting analyses for the 3-D motions of maser features. The 
procedure was performed on the basis of the least-square 
fitting of the observed kinematics to a radially-expanding 
flow model and in almost the same way as that applied 
to the W3 IRS 5 water masers (Paper I), which is not re- 
peated here. One difference is only the assumed speed of 
the radial expansion of a maser feature i, V(i), as a func- 
tion of distance of a maser feature from the originating 
point of the outflow, r^, which is expressed more simply 
as Vexp(«) = Vo(ri/ro) a , and where Vq is an expansion 
velocity at a unit distance ro, a is a power-law index in- 
dicating the apparent acceleration of the flow. We made 
the fitting step-by-step, excluding maser features having 
unreliably large positive or negative expansion velocities 



or distances from the outflow origin. After such several 
iterations, we used 68 proper motion data and obtained 
best solutions that are given in Table | 

Figure shows an estimated position of the outflow 
origin in the maser motions (plus), which coincides 
with the location of the SiO maser emission (filled 
square, EGHMM) within the position uncertainly 
(~100 mas). A systemic line-of-sight velocity of the 
flow is almost equal to that of the ambient molecular 
cloud (~56 km s _1 , e.g., S81; 000 [cite]cite.cox87Cox 
et al. (1987); 000 [cite]cite.rud90Rudolph et al. 
(1990); |00 [cite]cite.zha95Zhang & Ho(1995); 000 
[cite]cite.zha98Zhang, Ho, & Ohashi(1998); 000 
[cite]citc.oku01Okumura et al. (2001)) and roughly coin- 
cident with that of the SiO maser emission (~47 km s _1 ) 
within a velocity width of the cloud (~26 km s^ 1 ). 

Figure presents the maser feature motions projected 
onto three different planes. The maser kinematics indicate 
no rotation of the expanding flow, suggesting that ballistic 
motions predominate the kinematics. The best fit model 
and an expansion velocity plot against distance from the 
outflow origin (Figure 0) indicate that the expanding flow 
decelerates in the water maser region (r = 200-500 mas or 
1200-3000 AU from the outflow origin), where the expan- 
sion velocity decreases from ~90 km s _1 to ~50 km s _1 . 
This is a controversial case against these expanding flows 
that apparently exhibit the accelerations in water maser 
kinematics (Orion KL, 000 [cite]cite.gen81bGenzel et al. 
(1981)b; W 49N, §00 [cite]cite.gwi92Gwinn, Moran, & 
Reid(1992); W3 IRS 5, Paper I). 

3.3. W51 Main 

Figure 0a presents the angular distribution of water 
masers in W51 Main, which also revives those found by 
previous observations (000 [cite]cite.gen78Genzel et al. 
(1978); 000 [cite]cite.gen79Genzcl et al. (1979); G81) 
shown in Figure 0b. Four maser sub-clusters have 
been identified by G81: "Double Knot", "Middle High 
Velocity Cluster" , " Northern High Velocity Cluster" , and 
" Southern High Velocity Cluster" , all of which seem to be 
stable for at least 20 years. 

LLD identified a "cocoon" in the Double Knot, which 
is more clearly seen by superposing three maps of G81, 
LLD, and the present work around the coordinate (45, 
25) in unit of mas in Figure 0b. Assuming a distance to 
W51M of 6 kpc, this cocoon has an inner and an outer 
radii of approximately 12 AU and 60 AU, respectively. A 
rotation of the cocoon has been also identified by LLD, but 
it was not found in the present work, probably because of 
too small a number of maser features detected around the 
cocoon. 

On the other hand, the rcmaing maser sub-clusters have 
large position offsets up to 20 mas (120 AU at a distance 
of 6 kpc) between the previous and the present maps. 
Likely these are not true motions of sub-clusters but the 
" Christmas tree" effect due to appearance and disappear- 
ance of maser features during 20 years. Unlike the bow 
shock pattern seen in W 51N, no clear feature alignment 
was found except for the cocoon mentioned above. 
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Figure fj] presents the detailed 3-D motions of water 
masers around the four maser sub-clusters. Most of all 
mascr features with measured proper motions are red- 
shifted with respect to the systemic velocity of this re- 
gion (50-60 km s"\ e.g., G81; ©00 [cite]cite.cox87Cox 
et al. (1987); gOO [cite]cite.rud90Rudolph et al. (1990); ©00 
[cite]cite.zha95Zhang & Ho(1995); §00 [cite]cite.ho96Ho 
& Young(1996); §00 [cite]cite.zha98Zhang, Ho, & 
Ohashi(1998);§00 [cite] cite. okuO 1 Okumura et al. (2001)). 
Looking at the whole Doppler-velocity range (158 km s _1 
< Vlsr < — 58 km s^ 1 ). the kinematics of water masers 
is apparently random. On the other hand, LLD mea- 
sured 26 proper motions of water masers with a one-month 
time baseline and found a stream in the SW direction 
from the cocoon in a limited range of the Doppler veloc- 
ity (54 km s" 1 < Vlsr < 68 km s" 1 ). We note that the 
Southern High Velocity Cluster seems to independently 
have an expanding flow. Maser features on the north- 
east side of this sub-cluster (around the coordinate [40, 

— 100] in unit of mas in Figure 0) have proper motions to- 
ward the cocoon with velocities of 40-80 km s , probably 
which are not a flow contracting toward the cocoon but an 
expanding flow from the point around the coordinate (20, 

— 140) in unit of mas in Figure ^. The systemic Doppler 
velocity of the candidate originating this expanding flow 
is expected to be around Vlsr km s _1 on the ba- 
sis of the mean Doppler velocity of maser features in the 
sub-cluster. 

34. W51 South (e4) 

Figure § presents the 3-D motions of water maser 
features in W51 South, close to the 3.6-cm continuum 
source W51 c4 (§30 [cite]cite.gau93Gaume, Johnston, & 
Wilson(1993)) and 2-mm continuum source W 51 e8 (y00 
[citc]cite.zha98Zhang, Ho, & Ohashi(1998)). The distri- 
bution of water masers in this region seems to have been 
roughly stable and aligned in the north-west to south- 
east direction (G81). We found that the water masers 
exhibit marginally a Doppler-velocity gradient along the 
maser alignment and a systematic separation motion be- 
tween the red-shifted and the blue-shifted masers with 
bipolarity. We attempted a model fitting for the water 
masers using only a model such as that applied as Step 1 
for W51 North water masers, in which we esti mated only 
the originating point of the outflow (see Sect. 3.2.2 and 
Paper I). Table | gives its solution. 

4. Discussion 

4-1. Distance to the W 51A region 

The distance to the W51A region (G 49.5-0.4) has been 
adopted to be ~ 7 kpc on the basis of the statistical paral- 
lax of water-maser proper motions measured by G81 (c.f. 
S81). Because of an increase in measured proper motions 
and of reliable fitting of the maser kinematics to an ex- 
panding flow model in the present paper, a distance of ~6 
kpc should be adopted. 

As described in Sect. 3. 2.2] and Table §, we obtained a 
distance value of 6.1±1.3 kpc from a model fitting for the 



kinematics of W51N. The uncertainty permits a distance 
value of 7 kpc. When adopting 7 kpc, however, the model 
fitting puts the originating source of the W51N outflow 
at a point with a NW offset of ~0".l from that when 
adopting 6 kpc. As seen in Figures § and §, this offset is 
inconsistent with the suggestion that the originating point 
should coincide to the location of SiO masers (EGHMM). 
Moreover, although both of the NW and SE sub-clusters 
exhibit deceleration, unlike a simple deceleration as shown 
in Figure they arc not aligned on a simple one. 

We also attempted the statistical parallax method and 
obtained a distance value of 7.2±1.7 kpc but when includ- 
ing all measured proper motions in the W 51 North and 
West regions. When using only proper motions in the two 
dominant sub-clusters of W51N, an unreliably small value 
(~4 kpc) was obtained. We obtained a distance value of 
8.5±3.3 kpc from the kinematics of W51M. As described 
in Sect. 4.2 and Figures §and|4|, the 3-D maser kinemat- 
ics are heavily biased. Therefore, it is difficult to obtain a 
reliable distance value using the statistical parallax. Here 
we adopt the distance obtained by the model fitting. 

On the other hand, the distance to G49.5— 0.4 has been 
estimated to be 5.5 kpc on the basis of a 'far' kinematical 
distance adopting the Galactic constants: R®= 8.5 kpc 
and 9®= 220 km s _1 (c.f., [j]00 [cite]cite.cra78Crampton, 
Georgelin, & Georgelin(1978); §D0 [cite]cite.dow80Downes 
et al. (1980)), which is quite consistent with that we have 
obtained from the model fitting. Note that, for some of 
the Galactic water maser sources, the distances obtained 
from their water maser kinematics coincide quite well with 
their kinematical distances (e.g., W3 IRS 5, Paper I). 

4-2. Objective analyses of the maser kinematics 

In order to model-independently deal with the kinemat- 
ics of water masers in W 51 North and Main, we performed 
diagonalization for the velocity variance-covariance matri- 
ces obtained from velocity vectors of maser features (c.f., 
§00 [cite]cite.blo00Bloemhof (2000)). Table @ gives the 
obtained eigenvalues (velocity dispersions) and their cor- 
responding eigenvectors after the diagonalization. For ex- 
ample, for a bipolar outflow, an eigenvector correspond- 
ing to the largest eigenvalue indicates a major axis of 
the flow and others have two almost equal eigenvalues 
much smaller than the largest one. In the cases of 
both W51 North and Main, the obtained three eigen- 
values differ from each other by factors larger than two; 
the maser kinematics arc not isotropic but tri-axially 
asymmetric. Although the magnetic drag of bipolar 
flows might be applicable to explain such asymmetry 
(§00 [citc]citc.blo00Blocmhof (2000)), geometrical effects 
should be taken into account. 

A major eigenvector with the largest corresponding 
eigenvalue of the W51N kinematics (an inclination of 
16°±ll°and a position angle of — 41°±68°) is almost par- 
allel to the major axis of the W51N outflow, which can 
be confirmed in Figures | and |. EGHMM also obtained 
similar angles of the outflow axis on the basis of their 
observations of the W51N SiO masers. Having two eigen- 
values much larger than the one for the W51N kinematics 
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can be explained by the outflow having a large opening 
angle, which is much larger than that estimated from SiO 
masers (~25°, EGHMM) and by the flow being strongly 
blocked at some points. 

A geometrical effect such as that expected in W51N is 
also expected in the W51M kinematics. This can explain 
why most water masers are red-shifted with respect to 
the systemic velocity (Vlsr = 50-60 km s -1 , see Sect. 
|3.3| ) when supposing that the driving source of the water 
masers is located behind a dense molecular cloud core that 
blocks the flow material approaching us. This hypoth- 
esis is supported by the fact that the eigenvector with 
the smallest eigenvalue is almost along the line-of-sight. 
Note that possible multiple outflows along the sky plane 
as mentioned in Sect. 3.3 can also create such tri-axially 



asymmetric kinematics. 

4-3. Origins of the maser kinematics 

Water masers in massive-star forming regions are as- 
sociated with outflows from newly-formed massive stars. 
Therefore, we can estimate the locations of young massive 
stars as originating points of the outflows and, in some 
cases, elucidate the evolutionary statuses of the young 
stars and kinematical structures of the outflows them- 
selves as implied in Sect. 0. 

For the W51N flow, EGHMM proposed that the SiO 
maser should be associated with a conical flow that ro- 
tates (counter clockwise) about the axis of the expan- 
sion motion. On the other hand, we did not find such 
a helical motion in th e H2O maser kinematics as men- 
tioned in Sect. 



3.2.2 



Although H2O and SiO maser 
flows are driven from a common young star, the asso- 
ciating flows are likely streaming independently. Such 
a hybrid system was also proposed for the Orion-KL 
outflows (000 [cite]cite.gre98Greenhill et al. (1998)) and, 
therefore, might be a common characteristic of outflows 
from massive young stars. On the other hand, the W51N 
maser kinematics are heavily biased due t o po ssible block- 
ing of the flow as mentioned in Sect. 4.2. Adjacent 



young massive stars and their embedding dense cloud 
cores are expected as candidates of such blocking ob- 
jects because the W51N region is composed of several 
massive young stars within a very small volume (< 0.1 
pc) (e.g., 000 [cite] cite. wat98Watarai et al. (1998); 000 
[cite]cite.oka01Okamoto et al. (2001)). They may also 
create observed fast (<100 km s -1 ) random motions in 
the maser kinematics in fields far from the flow origin 
and between the cloud cores. It is difficult to create such 
random motions by means of shocks in the interstellar 
medium, which are triggered by adjacent HII regions, be- 
cause the HII regions themselves will expand much more 
slowly (Vexp < 30 km s _1 ). 

The existence of several massive stars in the W 51M 
is expected as mentioned in Sect. |3.3| . In fact, although 
a duration of H2O maser activity is quite short (<10 5 
yr, 000 [cite]cite.gen77Genzel & Downes(1977)), multi- 
ple driving sources in single water maser clusters have 
been found in some massive-star forming regions (e.g. W 3 
IRS 5, Paper I; Cep A, 000 [cite]cite.tor01Torrelles et al. 



(2001)). These imply that a mass infall rate in a star- 
forming region exhibiting such multiple outflows should 
be much higher than 10 -3 M®yr -1 . This lower limit of 
the rate is calculated using an assumed total mass of the 
star cluster larger than 100 M ® (in the case of W51N, 000 
[cite]cite.oka01Okamoto et al. (2001)) divided by the time 
scale within which this mass infalls before H2O masers are 
quenched by photoionization of H2O molecules due to the 
oldest massive star in the star cluster. 

4-4- Relative 3-D motions of water maser clusters in 
W51A 

In principle, it is possible to measure relative and sys- 
temic bulk motions between clusters of water masers in 
the same manner as the "in-beam astrometry" because 
they were observed simultaneously in a single beam of 
each telescope. In the present paper, the coordinates of 
all maser maps were precisely fixed with respect to the 
position-reference maser feature in W51N (W51N: 12002 
4). Table gives measured relative bulk motions. 

In practice, it is difficult to precisely estimate such bulk 
motions of the maser clusters because these motions are 
heavily disordered by random motions. Even for W51N, 
whose maser kinematics are relatively well known, the 
uncertainty of its systemic bulk motion is quite large 
(^40 km s _1 , sec Table 0) compared with the typical ve- 
locity width of a molecular cloud (< 10 km s _1 ). For 
W51M, only the mean proper motion of the W51M kine- 
matics is observable and adopted as a systemic bulk mo- 
tion of W51M; of course its uncertainty is also quite huge. 

Adopting the systemic bulk motions of W51 North 
and South obtained from model fittings (Table and 0) 
and the mean maser motion of W51 Main as the sys- 
temic bulk motion, we found a bulk motion separating 
W51 North and Main/South by -130 km s" 1 (c.f. 000 
[cite] cite. imaOHmai et al. (2001)). At first, the separation 
motion driven by pressures by the adjacent HII regions 
is considered. In fact, there are bright HII regions that 
are located 0.9 pc west from W51 Main/South (W 51 
IRS 1) and 0.1 pc north-cast from W51 North (W 51 
IRS 1/d) (e.g. 000 [cite]cite.gen82Genzel et al. (1982); 000 
[cite]cite.gol94Goldader & Wynn-Williams(1994)). On 
the other hand, the mechanical luminosity required for the 
bulk motions is calculated assuming a spherically expand- 
ing flow sweeping molecular clouds (maser regions) that 
have a density of molecular hydrogen, ?ih 2 ~ 10 6 enr 3 
(e.g., 000 [cite]cite.eli92Elitzur(1992)), in an expansion ve- 
locity to create the bulk motion, t>exp — 130/2 km s _1 , at 
a distance of r ~ 0.9 pc. However, the obtained luminos- 
ity, > 10 7 L® , is too huge to be produced only by stellar 
clusters in the W 51A region. The biases in the obtained 
systemic bulk velocities of the clusters are a ppar ently most 
likely to be generated as discussed in Sect. 4.2. Some ar- 
tificial effects on the coordinate drifts between the maser 
clusters are also expected: for example, unknown large 
offsets of telescope positions due to lack of geodetic VLBI 
observations with J-Net will result in the failure to resolve 
27r-n ambiguity and to connect fringe-phases through each 
of the observations. 
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5. Summary 

We have obtained data on the 3-D kinematics of three 
clusters of water masers in W 51A with VLBI monitoring 
observations composed of five epochs with an 8-month 
time baseline. The main conclusions of this paper are 
as follows. 

1. In W51N, a bipolar outflow is clearly exhibited by 
the maser kinematics, modelled using a radially expand- 
ing flow model. The position of the driving source of the 
outflow coincides with that of the SiO maser source within 
0".l. Expansion velocity is consistent with that obtained 
by EGHMM on the basis of the change in maser distri- 
bution for 15 years but decreasing with distance from the 
driving source from 90 km s _1 to 50 km s _1 . 

2. In W51M, random motions predominate the maser 
kinematics. Multiple outflows are considered. Unlike 
W51N, no systematic change in the maser distribution 
has been seen during 20 years due to the " Christmas tree 
effect" for maser appearance and disappearance. 

3. In W51S, a bipolar outflow was found marginally. 
The driving source of the outflow has a large offset (> 1") 
from both of the continuum sources W51 e4 and e8. 

4. The maser kinematics in W 51A are heavily biased 
due to possible concentration of young massive stars and 
their embedding clouds that destroy symmetry or system- 
atic motion of the outflows. Duration of active massive- 
star formation in each of the maser clusters is shorter than 
10 5 yr in W51A, which is suggested by the fact that the 
H2O masers associated with the different driving sources 
of outflows are simultaneously observed. 

5. The distance to W51A was estimated to be 6.1±1.3 
kpc on the basis of the model fitting method applied to the 
maser kinematics of W51N. On the estimated distance, 
the kinematical model for the W51N outflow is consistent 
with that previously proposed on the basis of VLA obser- 
vations (EGHMM) and exhibits simple deceleration in the 
flow. This distance value is smaller than that previously 
adopted (~7 kpc) but consistent with the kinematical dis- 
tance (~5.5 kpc). 

6. Measurements of relative 3-D bulk motions between 
the maser clusters have been attempted and showed an 
apparent separation motion between W 51 North and 
Main/South. To elucidate the true motions, random mo- 
tions and biases in the maser kinematics should be care- 
fully taken into account. Reliable kinematical models are 
indispensable to find accurately systemic bulk motions of 
maser clusters. 
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JU-24S10OO JD-245100* 

Fig. 1. Observed relative proper motions of water maser fea- 
tures in W51 North. Only the maser features detected at 
four or all five epochs are presented here. A number added 
for each proper motion shows the assigned one after the des- 
ignated name form "W51N:I2002" . A solid line indicates a 
lcast-square-fittcd line assuming a constant velocity motion. 
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Fig. 2. Dopplcr velocity distribution of water masers in the W51 North and West regions. A filled rectangle in the lower panel 
indicates the location of Silicon monoxide maser emission in W51N observed by EGHMM, which was estimated in the procedure 
mentioned in the main text. The distribution of the SiO masers obtained by EGHMM is also shown in the sub-frame. The opening 
angle of the SiO maser distribution has been extrapolated to the scale of the H2O maser emission (solid lines), which is the same as 
that of EGHMM. 
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Fig. 3. 3-D motions of water maser features in W51 North. A position, opening direction, and length of a cone indicate the position 
of a maser feature and the direction and amplitude of its proper motion, respectively. An inclination of a maser motion is calculated 
with respect to the systemic Doppler velocity of the parent molecular cloud (~56 km s — 1 ). Red and blue cones indicate receding 
and approaching features, respectively. An origin of the map is located at the position-reference maser feature. A plus at (—230 
mas, 130 mas) indicates an estimated position of the expanding flow in W51 North. A filled rectangle indicates the same as that in 
Figure 2. 
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Fig. 4. Estimated three-dimensional positions and motions of water maser features in W51 North. The positions and motions are 
with respect to those of the outflow origin determined by a model fit. The position of an arrow indicates that of a maser feature. 
Direction and length of an arrow indicates the direction and the magnitude of the maser motion, respectively. Top: Front view 
(XY-plane) of the positions and motions. Bottom left: Top view (XZ-plane). Bottom right: East-side view (Zy-plane). 
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Angular distance from the outflow origin (mas) 



Fig. 5. Radial expansion velocity of a water maser feature 
against a radial distance from the outflow origin in W51 
North after a model fit. According to the best-fit model, 
the expansion velocity is Vexp— 94 km s — 1 at r = 0".25 and 
Vexp^55 km s" 1 at r = 0".45. 
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Fig. 6. Same as Fig. 2 but in the W51 Main region, (a) Water masers found in the present work, (b) Maser distribution obtained 
by superposition of maps of G81 (observed in 1977-1979, red filled circle), LLD (observed in 1994, greed filled circle) and the present 
work (observed in 1999, blue filled circle). The coordinate of G81 was adopted. The double knot at (45, 25) in unit of mas was used 
as reference of the superposition. 
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Fig. 7. Same as Fig. 3 but for W51 Main. An inclination 
of a maser motion is calculated with respect to the systemic 
Doppler velocity of the parent molecular cloud (~55 km s ). 




—I i ■ i i I i i i i I i i i . I i i i ■ |_ 

SQ0 -500 -I0OG -15UU 



K.A. offset (mas) 

Fig. 8. Same as Fig. 3 but for W51 South. A plus at (230 
mas, —10 mas) indicates an estimated position of the expand- 
ing flow in W51 South. 
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Table 1. Status of the telescopes, data reduction, and resulting performances in the individual epochs of the J-Net observations. 





Epoch in 






Reference 


l-cr level 


Synthesized 








Observation 


the year 


Duration 


Used 


velocity 2 


noise 


beam 3 


Features 




code 


1999 


Or) 


telescopes 1 


( km s _1 ) 


(Jy beam - ) 


(mas) 


N 


M 


S 


j99052 .... 


February 21 


11 


M, b, k 


37.74 


9.3 


7.4x2.4, -38.4° 


137 


51 


22 


j99092 .... 


April 2 


10 


M, b, n 5 , k 


34.67 


3.8 


8.8x2.9, -51.6° 


100 


50 


17 


j99126 .... 


May 6 


10 


M, b, n, k 


35.72 


2.9 


8.1x3.0, -41.3° 


171 


88 


23 


j99157 .... 


June 6 


10 


M 6 , b 6 ' 7 , n 6 , k 6 


35.51 


4.0 


7.2x2.6, -54.2° 


61 


33 


12 


j99294 .... 


October 21 


11 


M, b, k 


53.95 


22.0 


7.6x2.4, -38.2° 


47 


25 


6 



1 Attending telescopes; M: the 10-m telescope at Mizusawa, b: the 34-m telescope at Kashima, n: the 45-m telescope 
at Nobeyama, k: the 6-m telescope at Kagoshima. 

2 Velocity channel used for the phase reference in data reduction. 

3 The synthesized beam made in natural weight; major and minor axis lengths and position angle. 

4 Number of detected maser features in the W51A region (North, Main, and South [c4]). 

5 Ceasing operation for 7 hr due to strong winds. 

6 High system temperature due to bad weather conditions. 

7 Data heavily affected by unstable local-frequency signals. 
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Table 2. Parameters of the water maser features identified by proper motion toward W51 North. 



Offset 1 * Proper motion' 3 Radial motion 4 Peak intensity' 



Maser Feature 1 




(mas) 








(mas 


_ 1 

yr 


) 






(km s" 


- 1 ) 








Jy beam 1 


) 






(W51N:I2002) 


R.A. 


decl. 




/'.i 






/'.« 


a fly 


V z 




j99052 


j990092 


j99126 


j99156 


j99294 


1 


79 


47 


-22 


41 


1 


00 





72 





59 





76 


-60.28 


2 


85 




19 


97 


46 


08 








2 


88 


04 


47 


52 


1 


92 





73 


1 


29 





51 


-54.49 


1 


33 




12 


82 


37 


50 


7 


03 




3 


89 


08 


47 


20 





84 





29 


1 


24 





28 


-52.22 


2 


00 


892.23 


130 


49 


211 


22 


8 


83 


156.40 


4 





00 





00 





00 





36 





00 





35 


-47.71 


1 


85 


110.23 


8 


71 


9 


17 


2 


05 


144.14 


5 


72 


00 


-0 


00 


1 


00 


1 


10 


-1 


47 


1 


32 


-46.45 


1 


80 




5 


99 


12 


98 








6 


72 


39 


-1 


23 





38 


1 


82 


-1 


38 


2 


37 


-44.81 


2 


74 


100.20 


5 


41 












7 


11!) 


36 


31 


27 





00 


1 


81 


1 


52 


1 


11 


-39.86 


1 


(iO 




8 


28 


10 


32 








8 


129 


20 


7 


02 


1 


00 





46 





80 





92 


-34.65 


2 


25 


45.38 


5 


57 


80 


25 








9 


125 


34 


31 


14 





08 


2 


38 





76 


2 


45 


-34.24 


1 


79 




8 


50 


5 


00 








10 


124 


67 


38 


42 





50 


1 


81 


1 


20 


1 


8,3 


-32.96 


1 


80 


190.68 


5 


24 












11 


125 


02 


38 


57 


-0 


70 





08 


2 


14 





00 


-31.83 


1 


54 




2 


80 






5 


58 


30.60 


12 


125 


45 


30 


24 


-4 


20 


1 


10 


1 


05 


1 


60 


-31.41 





74 








3 


03 






14.90 


13 


129 


38 


8 


30 





55 


2 


05 


-0 


08 


2 


42 


-31.31 


1 


47 




2 


73 


13 


37 


2 


83 




14 


125 


69 


38 


30 


-0 


38 


1 


05 


2 


72 


1 


27 


-30.37 


1 


54 




1 


71 


17 


69 


8, 


83 




15 





63 


-7 


08 


-1 


03 





50 


-0 


02 





48 


-28.94 


1 


37 


89.45 


6 


00 


20 


87 


2 


48 




16 


79 


72 


10 


04 


1 


73 


1 


41 


1 


40 


1 


50 


-28.66 


1 


00 




!) 


01 


3 


00 








17 


82 


41 


-25 


44 


-4 


08 


2 


09 


3 


00 


2 


10 


-26.16 





84 




1 


98 


4 


03 








18 


130 


33 


6 


28 





07 





50 


1 


36 





07 


-25.81 


1 


80 




2 


77 


12 


12 


3 


95 


25.10 


19 


43 


26 


-17 


39 





71 


1 


56 





15 


1 


33 


-25.33 





04 


44.58 





17 












20 


126 


16 


4 


00 


-0 


12 





21 





88 





29 


-24.58 


1 


10 


61.40 






15 


29 


2 


48 


46.81 


21 


66 


70 


78 


80 





24 


1 


81 


1 


06 





90 


-24.21 


1 


37 


247.27 


10 


80 












22 


66 


88 


78 


89 


-0 


on 





47 





02 





38 


-23.14 


1 


89 


151.73 


29 


84 


00 


22 





21 




23 


109 


27 


-25 


40 


-2 


10 


2 


14 


-1 


49 


2 


44 


-20.22 


1 


37 


188.58 


2 


48 












24 


125 


98 





70 


-0 


on 





29 





08 





27 


-20.14 


2 


15 


64.66 


10 


81 


01 


50 


14 


81 


46.37 


25 


-19 


76 


-53 


79 


-1 


91 


1 


20 


-0 


05 


1 


05 


-20.10 





04 


155.72 


7 


09 












26 


109 


32 


-25 


55 


-1 


13 





25 


-1 


30 





60 


-16.36 


2 


82 


649.10 


29 


12 


70 


13 








27 


21 


90 


-66 


78 


-1 


54 





71 


-1 


50 





78 


-15.71 


1 


26 


96.70 


18 


81 


37 


00 








28 


128 


27 


-0 


23 


-0 


70 





20 


1 


64 





35 


-14.45 


2 


11 








139 


10 






37.50 


29 


398 


00 


-628 


58 


-1 


94 





34 





02 





73 


-14.43 


2 


85 


457.59 


40 


20 


04 


05 


15 


05 




30 


109 


33 


-25 


85 


-0 


30 





57 


-0 


65 





47 


-14.34 


2 


10 


427.39 


45 


80 


7 


99 








31 


104 


82 


-27 


55 





33 





66 


-1 


78 





34 


-12.80 


3 


74 


339.00 


37 


72 


108 


41 


30 


10 




32 


397 


93 


-628 


08 


-3 


37 


2 


23 


3 


41 
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-413 
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Water Masers in W51A 



Table 2. (Continued) 



Mascr Fcatu 



(W51N:I2002) 



Offset 2 
(mas) 



Proper motion 
(mas yr~ 1 ) 



Radial motion 
(km s - 1 ) 



Peak intensity 
(Jy beam -1 ) 



R.A. 



dccl. 



V z 



j99052 



76 
77 

78 

79 

80 

81 

82 

83 

84 

85 

86 

87 

88 

89 

90 

91 

92 

93 

94 

95 

96 

97 

98 

99 

100 

101 

102 

103 

104 

105 

106 

107 

108 

109 

110 

ill 

112 

113 

114 

115 

110 

117 

118 

11!) 

120 

121 

122 

123 



-401.66 
835.32 
-427.69 
857.86 
-367.68 
-388.12 
-401.82 
-423.16 
-17.07 
-426.36 
-428.42 
-402.29 
-401.42 
-431.38 
-389.02 
-394.69 
-427.22 
-413.70 
-393.76 
-402.97 
-435.49 
-402.62 
-422.76 
-417.88 
-414.21 
-422.90 
-414.44 
-424.43 
-422.89 
-424.90 
-76.82 
-419.15 
-425.23 
-428.83 
-428.82 
-370.82 
-409.02 
-423.54 
-373.63 
-421.19 
-386.01 
-420.92 
-383.90 
-375.03 
-425.92 
-420.55 
-422.97 
-421.18 



321.88 
680.28 
358.49 
619.66 
375.19 
384.99 
405.92 
325.87 
136.14 
346.63 
364.67 
324.52 
327.64 
365.57 
384.48 
385.29 
365.44 
279.51 
384.19 
324.56 
358.26 
323.02 
231.79 
358.94 
280.46 
231.69 
281.71 
226.00 
231.03 
203.86 
-40.30 
224.97 
263.61 
214.15 
215.68 
229.75 
230.41 
214.87 
228.51 
222.93 
236.21 
222.06 
232.24 
227.21 
215.03 
223.46 
216.81 
214.93 



-3.83 
-3.19 
-2.99 
-0.45 
-3.60 
-1.56 
-1.37 
-1.34 
-1.63 
-3.34 
-3.59 

0.59 
-3.07 

1.49 
-3.34 

3.36 
-4.36 
-5.31 

1.89 
-1.13 
-3.51 
-0.31 
-2.72 
-3.62 
-3.39 
-3.20 
-2.51 
-4.09 
-4.37 
-4.95 
-2.14 
-4.49 
-3.21 
-3.85 
-3.92 
-4.01 
-3.60 
-1.41 

0.40 
-3.91 
-4.54 
-5.31 
-3.19 
-5.16 
-2.17 
-3.02 
-4.24 
-4.60 



0.56 
0.35 
0.92 
2.83 
0.99 
0.87 
1.39 
1.12 
1.16 
0.89 
1.42 
2.04 
3.12 
2.10 
2.43 
2.71 
3.05 
3.34 
2.75 
1.66 
0.93 
3.23 
1.33 
0.83 
1.00 
0.51 
0.81 
3.29 
1.07 
1.07 
0.80 
2.59 
0.19 
0.76 
0.26 
2.72 
0.96 
3.42 
3.59 
0.58 
0.47 
1.66 
3.54 
2.32 
0.40 
0.30 
0.87 
0.83 



0.85 
3.17 
3.05 
1.33 
4.34 
1.58 
2.46 
0.49 
2.19 
4.16 
4.04 
1.19 
4.04 
2.26 
1.73 

-0.71 
3.12 
9.68 

-0.94 
0.83 
3.09 
0.99 
1.35 
3.89 
5.44 
1.32 
3.70 
0.85 
3.81 
1.76 
0.11 
3.21 
3.97 
9.80 
1.16 
2.79 
2.22 
0.92 
8.19 
3.24 
4.13 
3.20 
3.72 
2.25 
0.67 
1.51 
2.96 
0.86 



0.76 
0.29 
1.06 
2.69 
1.09 
1.15 
1.24 
1.21 
1.01 
0.84 
1.20 
1.33 
3.44 
1.64 
1.56 
3.22 
2.33 
3.34 
3.55 
2.78 
0.99 
3.47 
1.52 
0.84 
1.33 
0.78 
1.16 
4.27 
1.61 
1.10 
0.78 
2.99 
0.29 
1.29 
0.61 
3.76 
1.03 
3.09 
2.85 
0.65 
0.66 
1.56 
4.52 
2.07 
0.48 
0.62 
1.29 
1.05 



14.12 
14.40 
14.79 
14.90 
15.21 
15.27 
15.66 
16.14 
16.43 
17.57 
19.87 
19.93 
20.02 
20.07 
20.72 
21.83 
22.12 
22.13 
22.99 
23.91 
24.19 
24.48 
24.85 
25.29 
26.73 
27.35 
27.42 
29.03 
30.16 
32.26 
33.95 
35.42 
36.75 
41.94 
44.43 
49.11 
52.89 
54.79 
55.71 
55.79 
57.96 
58.81 
59.37 
61.86 
62.05 
64.05 
67.83 
69.41 



1.21 
2.42 
0.98 
1.58 
0.84 
0.84 
0.94 
1.78 
1.89 
1.12 
1.15 
1.47 
1.47 
1.68 
1.57 
1.15 
1.05 
0.94 
1.15 
2.74 
1.37 
1.57 
1.90 
1.55 
0.98 
1.68 
0.67 
0.84 
1.15 
0.63 
1.40 
2.32 
3.21 
1.90 
2.32 
1.15 
1.58 
0.37 



.89 
.47 
.12 
.75 
.17 
.06 
1.05 
3.72 
2.11 
2.95 



37.37 
140.49 
18.33 

26.70 
37.70 
33.46 
88.10 

141.38 
37.96 
18.83 

189.65 

62.64 
44.34 

39.54 

75.53 

42.87 
22.56 
49.16 
60.51 
26.63 
82.30 
15.40 

22.10 
18.90 
25.26 



28.00 
256.34 



10.90 
26.49 
15.24 



1701.62 
18.70 
124.73 



j990092 
7.06 
6.36 
11.57 
2.22 

7.13 

15.43 
5.29 



7.00 

4.53 
5.52 

6.42 
7.02 
11.65 



6.97 
56.67 
10.95 
4.64 
3.96 



3.63 
3.32 
7.23 

62.45 



1.56 
2.69 



3.14 

2.88 
2.62 
27.69 
7.35 
9.66 



j99126 
10.14 
12.34 
11.98 
7.51 
7.01 
5.21 
5.25 



6.77 
6.00 

8.74 

29.47 
5.95 



j99156 j99294 
2.81 

... 49.54 



8.62 

2.11 
13.45 

5.76 
15.01 
16.69 
21.30 
11.55 

4.63 



3.22 



7.46 



61.21 11.47 
8.19 



10.34 
1.51 

15.73 
9.93 
3.33 
2.38 
1.98 
6.78 
9.05 

32.30 

17.86 
115.72 
3.01 

12.79 



2.63 
2.30 



4.01 
1.86 



2.72 

8.28 



7.40 
2.49 



2.13 
1.78 



36.98 
59.33 



1 Water mascr features detected toward W51 North. The feature is designated as W 51N:I2002 N, where N is 
ordinal source number given in this column (12002 stands for sources found by Imai et al. and listed in 2002). 

2 Relative value with respect to the location of the position-reference maser feature: W51N:I2002 4- 

3 Relative value with respect to the motion of the position-reference maser feature: W51N:I2002 4- 

4 



Relative value with respect to the assumed systemic velocity of V^gp^ 
5 Peak intensity at five epochs. 



56.0 km s 



18 
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Table 3. Same as table 2 but toward W51 Main. 



Offset^ Proper motion^ Radial motion 4 Peak intensity^ 



Maser Feature 1 
(W51M:I2002) 




(mas) 








(mas 


yr~ x ) 








(km 


s- 1 ) 










(Jy beam 






R.A. 


dccl. 




!>.r 






/ 1 u 




V 2 




j99052 


j990092 


j99126 


j99156 


J99294 


1 


1289 


70 
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9 
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1 


18 
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51 


2 


29 
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94 


38 
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2 
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3 
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71 


9 


TO 


2 


07 


8 


i u 
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1 
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4 


— 28 


80 
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18 
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77 


1 


43 


1 


11 
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2.33 


1 


37 
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56 
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1 


15 














174.50 


2312.00 
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08 
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66 


3 


97 


2 


44 


1 


72 


2 


40 


30.31 


1 


05 


18 


00 


7 


23 










29 


30 


02 


-104 


23 


1 


77 





30 


4 


34 





47 


32.82 


1 


16 


45 


13 


16 


31 






9.12 


202.00 


30 


-15 


07 


-66 


38 


2 


80 





35 


3 


12 





48 


32.89 


1 


10 










10 


41 




132.75 


31 


-24 


22 


-66 


13 


1 


54 





08 


2 


44 





94 


33.65 


1 


26 


35 


75 






32 


44 






32 


32 


54 


-106 


43 





75 


1 


41 


4 


15 


1 


20 


36.84 





91 


100 


00 


14 


81 






4.47 




33 


34 


57 


-106 


32 


1 


64 





35 


5 


70 





32 


38.88 


1 


70 


325 


00 


92 


50 






16.63 


303.00 


34 


-30 


00 


-167 


34 


-1 


27 


2 


51 


3 


28 


1 


70 


39.76 


1 


54 






36 


81 


28 


94 


7.16 




35 


38 


0!) 


-107 


99 


4 


20 


2 


37 


3 


95 


1 


52 


41.31 


1 


05 


86 


62 


40 


37 










36 


37 


25 


-108 


55 


2 


47 


3 


75 


4 


08 


9 


63 


42.74 


2 


42 


107 


00 


107 


50 










37 


37 


11 


-109 


39 





70 


1 


34 


4 


42 


1 


15 


43.63 


2 


42 














8.94 


1278.00 


38 


40 


77 


87 


07 





13 


1 


92 


4 


00 


1 


01 


46.48 


1 


.37 






22 


81 


37 


12 






39 


28 


23 


-115 


30 


2 


37 


2 


83 


2 


89 


2 


26 


50.73 





84 


25 


25 


4 


84 










40 


-26 


72 


-164 


19 


6 


17 


4 


00 


3 


11 


3 


71 


52.64 


1 


47 


16 


80 


7 


25 










41 


-27 


11 


-162 


25 





48 


1 


00 


4 


20 


1 


68 


56.36 


1 


97 


09 


50 


19 


09 


7 


72 






42 


39 


75 


40 


08 





47 


1 


00 


5 


50 


1 


34 


59.05 


1 


62 


22 


87 


12 


44 


12 


02 






43 


-28 


60 


-81 


84 


7 


03 


1 


07 





00 





08 


61.16 


1 


01 


106 


00 


22 


88 


12 


37 






44 


-2 


27 


-140 


11 


4 


19 


3 


73 


1 


43 


2 


32 


63.03 


1 


37 


110 


00 


4 


66 










45 


32 


00 


-91 


38 





87 





57 


3 


03 





71 


65.62 


1 


15 










4 


87 




209.00 


46 


14 


72 


-142 


01 


1 


17 





90 


1 


54 


1 


09 


67.79 


2 


00 


111 


00 


35 


75 


33 


94 


4.16 




47 


27 


03 


-370 


19 





80 


3 


47 


1 


30 


2 


21 


80.99 


1 


48 






11 


19 


52 


75 






48 


29 


61 


-377 


16 


-4 


22 


3 


30 


7 


13 


3 


76 


101.63 





84 






8 


53 


8 


56 







1 Water maser features detected toward W 51 Main. The feature is designated as W51M:I2002 N, where N is 
ordinal source number given in this column (12002 stands for sources found by Imai et al. and listed in 2002). 

2 Relative value with respect to the location of the position-reference maser feature: W51M:I2002 15. 

3 Relative value with respect to the motion of the position-reference maser feature: W51N:I2002 4- 

4 Relative value with respect to the assumed systemic velocity of Vlsr = 55.0 km s _1 . 

5 Peak intensity at five epochs. 
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Table 4. Same as table 2 but toward W51 South (e4). 



Offset** Proper motion' 5 Radial motion 4 Peak intensity^ 

Mascr Feature (mas) (mas yr~ ) (km s - ) (Jy beam - ) 



(W51S:I2002) R.A. decl. /x x cr^ x fi v ay. v V. AV- j99052 j990092 j99126 j99156 j99294 

1 281.54 -261.33 2.03 4.18 1.84 7.78 -106.14 1.05 ... ... 15.10 3.41 

2 259.42 -255.53 2.56 1.46 0.55 2.15 -91.88 1.89 ... ... 39.82 12.82 

3 -1049.26 -145.00 -5.00 1.85 0.35 2.41 -7.92 1.68 ... ... 383.00 6.41 

4 -29.80 -6.27 -1.82 0.28 2.43 0.29 -7.26 1.68 ... ... 352.75 51.22 825.00 

5 -2.93 1.73 -2.28 1.18 2.58 1.00 0.80 1.19 210.00 ... 60.53 6.88 

6 0.00 0.00 -1.63 0.55 2.38 0.33 2.01 1.51 275.00 139.69 108.37 16.27 390.50 

7 -323.67 37.11 -1.85 0.55 1.84 0.52 2.50 0.70 99.50 ... 28.38 5.67 

8 -315.17 36.83 -0.57 0.79 2.74 0.77 3.99 0.74 38.78 ... ... 5.01 

9 -269.53 17.89 -2.16 0.20 2.29 0.19 4.81 1.55 1016.00 ... 494.00 70.31 

10 -828.27 572.16 2.41 5.71 3.26 3.58 15.00 1.37 ... ... 38.56 6.78 



1 Water maser features detected toward W 51 South. The feature is designated as W51S:I2002 N, where N is the 
ordinal source number given in this column (12002 stands for sources found by Imai et al. and listed in 2002). 

2 Relative value with respect to the location of the position-reference maser feature: W51S:I2002 6. 

3 Relative value with respect to the motion of the position-reference maser feature: W51N:I2002 4- 

4 Relative value with respect to the assumed systemic velocity of V^gR ~ 59.0 km s _1 . 

5 Peak intensity at five epochs. 



Table 5. Best-fit models for the maser velocity field of the W51 North outflow 



Parameter 


Step l 1 


Step 2 Z 


Offsets 


Velocity: 

^ 0x 3 ( km s- 1 ) 


-52±18 


-50±19 


Vq/ ( km s" 1 ) 
Vq^ (kms- 1 ) 


22±12 
5 


14±34 
2±7 


Position: 






xq (arcsec) .... 


-0.28±0.06 


-0.23±0.08 


Uq (arcsec) .... 


0.11±0.05 


0.13±0.10 


Velocity field 


Radial outflow: 






V\ ( km s 1 arcsec Q ) ... 6 


27±6 


a 


6 


-0.9±0.2 


Distance d (kpc) . . 


6.0 7 


6.1±1.3 


RMS residual VS^ 


16.23 


3.67 



1 Assuming independent expansion velocity of maser features. 

2 Asumming a common expansion velocity as a function of distance from the outflow origin as shown in the main 
text. 

3 Relative value with respect to the position-reference maser feature. 

4 Relative value with respect to Vlsr = 56 km s _1 . 

5 Step 1 assumes the systemic radial velocity: Vq z = 0.0 km s _1 . 

6 The solution determines a radial outflow velocity Vexp (i) independently for each feature with a proper motion. 

7 Distance is completely covariant with the Z[ and Vexp(i) and cannot be determined: d=7.Q kpc. 
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Table 6. Best-fit model for the maser velocity field of the W51 South (e4) outflow 



Velocity offset: 



(kms- 1 ) 



Vi 



oy; 



( km s ) 



y nz 2 ( km s- 1 ) 



83±15 
71±8 
3 



Position offset: 



xq (arcsec) . . 
j/q (arcsec) . . 
Distance d (kpc) 



0.23±0.14 
-0.01±0.02 
6.0 4 



RMS residual VtP 



9.32 



1 Relative value with respect to the position-reference maser feature in W 51N (W51N: 12002 4)- 

2 Relative value with respect to Vlsr = 59 km s -1 . 

3 Assuming the systemic radial velocity: Vq z = 0.0 km s — 1 . 

4 Distance is completely covariant with the z\ and Vexp(i) and cannot be determined: d=6.Q kpc. 



Table 7. Eigenvalues and eigenvectors of velocity variance-covariance matrixes of the maser kinematics in W51 North and Main 
after diagonalization. 



Eigenvalue 




(km 2 s- 2 ) 


Normalized eigenvector 


W51N 4759 


(-0.6332, 0.7141, 0.2984) 


2084 


(0.6819, 0.6971, -0.2213) 


406 


(0.3660, -0.0634, 0.9284) 


W51M 7169 


(0.9214, -0.3703, -0.1182) 


2381 


(0.3887, 0.8778, 0.2800) 


653 


(0.0001, -0.3039, 0.9527) 



Table 8. Relative bulk motions of three maser clusters with respect to the position-reference feature in W 51 North. 



Vy W Used 

( km s -1 ) ( km s _1 ) ( km s _1 ) proper motions 



W51 North (average) 


-52±54 2 


36±58 2 


60±35 2 


92 


(model fit) 


-50±19 


14±34 


58±7 


68 


W51 Main (e2) (average) 


61±63 2 


91±45 2 


70±40 2 


37 


W51 South (e4) (average) 


-12±66 2 


61±25 2 


29±54 2 


10 


(model fit) 


71±15 


59±8 


59 3 


10 



1 Velocity with respect to the local standard of rest. 

2 Uncertainty was estimated from the standard deviation of motions. 

3 Assuming Vlsr = 59 km s _1 . 



